Rheology of human blood plasma: Viscoelastic versus Newtonian behavior 
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We investigate the rheological characteristics of human blood plasma in shear and elongational 
flows. While we can confirm a Newtonian behavior in shear flow within experimental resolution, 
we find a viscoelastic behavior of blood plasma in the pure extensional flow of a capillary break-up 
rheometer. The influence of the viscoelasticity of blood plasma on capillary blood flow is tested 
in a microfluidic device with a contraction-expansion geometry. Differential pressure measurements 
revealed that the plasma has a pronounced flow resistance compared to that of pure water. Supple- 
mentary measurements indicate that the viscoelasticity of the plasma might even lead to viscoelastic 
instabilities under certain conditions. Our findings show that the viscoelastic properties of plasma 
should not be ignored in future studies on blood flow. 

PACS numbers: PACS 83.50. Jf 83.60.Wc 87.19.U- 



Blood is a complex fluid that consists of a suspension 
of blood cells in a liquid plasma which contains mostly 
water as well as proteins, mineral ions, hormones, and 
glucose. In humans, red blood cells (RBC) are the most 
abundant type of cells in whole blood with a concentra- 
tion of approximately 45% by volume. Because of this 
high RBC concentration, it is often believed that rheo- 
logical behavior of whole blood is mostly determined by 
the presence of the RBCs. Blood exhibits shear thin- 
ning and at low shear rates the red blood cells form 
reversible aggregates (rouleaux) that are broken up at 
high shear rates [Ij. The rouleaux formation is caused 
by the plasma proteins (most likely due to a depletion 
effect) but the plasma solution, which is approximately 
92% water, is currently believed to be Newtonian [lH3]. 
For example, experiments by Copley and King [4] and 
Jaishankar, Sharma and McKinley [5j show that the non- 
Newtonian and viscoelastic effects found in plasma and 
BSA (Bovine Serum Albumin) solutions in shear flows 
can be attributed to the surface layer of plasma pro- 
teins present at the liquid-air interface; that is, these 
non-Newtonian effects are surface rather than bulk ef- 
fects. These surface effects can be suppressed by the 
addition of a small amount of surfactants to the protein 
solutions [5J. 

Much work has been devoted to the understanding of 
laminar blood flow as well as the appearance of flow insta- 
bilities that may appear under pathological conditions, 
e.g. near an aneurysm or blockage [6, 7J. In general, one 
has to distinguish between inertia and viscoelastic driven 
instabilities. An elegant way to separate these two kinds 
of instabilities is the observation of vortices in a con- 
traction flow device. Inertial effects lead to downstream 
vortices whereas upstream vortices can be observed for 
viscoelastic fluids ^ . This approach has been used in 
a recent work by Sousa et al. [2j which gives an overview 
of the flow of two common blood analog solutions in 



contraction-expansion microfluidic devices. The two in- 
vestigated plasma replacement solutions with the same 
shear but different elongational viscosity revealed a re- 
markably different flow behavior which shows that it is 
not sufficient to consider only the shear viscosity in or- 
der to build a plasma analog solution. A pronounced 
viscoelastic behavior in elongational flow has been found 
in several other bio-fluids, too, e.g. for saliva fiO\ [TT] or 
DNA solutions. 
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Figure 1. Snapshots of the capillary bridges in the CaBER 
experiment. Images are 833 |J,m x 1344 |j,m. 

In this manuscript, we investigate the rheology and 
flow behavior of human blood plasma in both shear and 
extensional flows. We find that blood plasma shows sig- 
nificant viscoelastic effects in elongational flow in both (i) 
a capillary breakup extensional rheometer (CaBER, Fig- 
ur^ and (ii) a microfluidic contraction-extension device. 
In addition a synthetic (polymeric) plasma replacement 
solution is identified that matches the shear and exten- 



2 




0.000 



0.002 p-<>-<>-^>-<>-<>-0-<>-<>-<>-<>-<>-<>-<>-<>-0 
0.001 



M < 2|jNnn 



10 



100 



1000 



Figure 2. (Color online) Shear rate dependent viscosity 77(7) 
data at T = 20° C obtained by the rotational shear measure- 
ments. The grey shaded area indicates the regime below the 
resolution of the rheometer. 
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r(°c) 


?7o(mPa s) 


Ac(ms) 


Whole Blood 


37 


16.9 


7.8 ±0.6 


Plasma 


20 


1.95 


2.6 ±0.2 


Plasma 


37 


1.34 


1.5 ±0.2 


P£;o500i5/85 


20 


2.47* 


15.7 ±0.2 


P£;O50i5/85 


20 


1.40 


2.6 ±0.1 


P£;O5055/45 


20 


7.88 


7.9 ±0.2 


P£;O50i5/85-i/ct50 


20 


8.03* 


1.84 ±0.04 


BSAH20 


20 


1.24 


Newt. 




37 


7.6* 


Newt. 


H2O 


20 


0.97 


Newt. 



Table 1. (Zero) shear viscosities 770 and CaBER relaxation 
times Ac of the sample solutions. Values marked with * are 
obtained by a fit based upon the Carreau model [12]. 



sional rheology of human plasma. 

Sample solutions: Blood plasma used in the experi- 
ments was is taken from three healthy donors and stored 
in test tubes using anticoagulants (EDTA or Citra). 
Plasma is obtained by centrifuging the blood samples 
at 4000 rotations per minute for 5 minutes immedi- 
ately after withdrawal. The separated transparent liquid 
phase is centrifuged again and finally the pure plasma 
phase is extracted. All measurements are performed 
on the day of donation. Since the main protein in the 
plasma is albumin, an aqueous solution of 45 mg/ml 
albumin (BSA, Polysciences) is prepared for compari- 
son. The plasma replacement solutions are prepared by 
adding small amounts of polyethylene oxide (PEO, Sigma 
Aldrich, = 4-10^ g/mol) to a Newtonian solvent. 
Three PEO solutions are prepared: (i) a 50 ppm PEO 
in 15/85wt% glycerol/ water solution (P£^O50i5/85), (ii) 
a 500 ppm PEO in 15/85wt% glycerol/ water solution 
(P£;O500i5/85), and (iii) a 50 ppm PEO in a 55/45^t% 
glycerol/water solution (P£^O5055/45). Two RBC sus- 
pension at 50vol% are also used, one in physiological 
buffer (PBS, Phosphate Buffered Saline, Invitrogen) with 
15wt% glycerol (PP5'i5/85-i^ct5o) and the other under 
the same condition but with the addition of 50 ppm of 
PEO (P£^050i5/85_i:/ct5o)- Finally, whole blood and dis- 
tilled water are used as reference fluids. 

Steady shear rheology is performed in a rheometer 
(Thermo Scientific, HAAKE MARS II) at shear rates 
ranging from 7 = 1 to 1000 in rate controlled 
(CR) mode. A double-cone DC60/1° geometry is used in 
order to avoid the effect of surface viscoelasticity on the 
measured shear viscosity r]{j). To reduce the instrumen- 
tal noise, a special protocol is applied. For each shear 
rate, the torque signal is integrated over at least a full 
revolution of the cone. The minimal resolvable torque is 
Tmin = 2 iLiNm which limits the reliable minimal shear 
rate. Above, the accuracy is A77 = ±0.1 mPas. The 



temperature is kept constant at T = 20 ± 1°C. 

Figure [2] shows the shear viscosity of blood plasma at 
room temperature {RT = 20° C), along with viscosity 
data of whole blood (average of three donors), two di- 
lute PEO solutions, and a reference curve of pure water. 
Within the reliable data range (above the grey shaded 
area), the plasma sample as well as the 50 ppm PEO 
in 15wt% glycerol- water solution show Newtonian be- 
havior, that is constant shear viscosity rj. By contrast 
the addition of a physiologic amount of RBCs (Hct50%) 
significantly increases the viscosity of the PEO solu- 
tion. The solution exhibits shear thinning and the flow 
curve can be fitted to the Carreau model [12]. The vis- 
cosity of full blood is higher than the viscosity of the 
P£^O50i5/85-//c^50 solution. It is known that plasma 
proteins lead to a depletion induced formation of aggre- 
gates of RBCs at low shear rates, and this effect might 
be less pronounced for the PEO solutions [Ij. Table [l| 
summarizes the estimated zero shear viscosities 770- We 
should mention that our determination of 770 for whole 
blood is of limited use only because blood is a yield stress 
fluid (yield point was found to be approximately 5 mPa) . 

The relaxation time of the fluids is obtained using a 
custom made capillary breakup extensional rheometer 
(CaBER) setup [HI [ll]. The thinning of the capillary 
bridge is recorded using a high speed camera and 4x mag- 
nification objective (Pentax, PLAN 4x/0.10, WD=30). 
We find that the solutions have relatively short relax- 
ation times, and a slow retraction method is used to al- 
low for a quantitative evaluation of the data |15j. For a 
non-Newtonian fluid, the initial inert io-capillar thinning 
of the capillary bridge is followed by a formation of a 
straight filament, which thins exponentially in time. The 
relation h{t) = ho ■ ex.-p { — {t — to)/ Xc) gives the char- 
acteristic CaBER relaxation time Ac, where h{t) is the 
diameter of the filament at the time t during the exponen- 
tial thinning process, which starts at time to and width 
ho respectively. According to the Oldroyd-B model, a 
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Figure 3. (Color online) Width h{t) of the capillary bridge 
as a function of time t — tfinai- Markers on the right vertical 
axis refer to the images shown in Figure [l] Straight lines are 
exponential fits. 



simple approach leads to an estimate for the polymer re- 
laxation time A = Ac/3 [16} flT]. 

Figure [T] shows a series of images of the thinning pro- 
cess of the capillary bridges of the plasma as well as the 
BSA and P£^O50i5/85 solutions. The capillary bridge of 
the BSA solution shows the same dynamics as the dis- 
tilled water (not shown), within experimental resolution. 
In this case, the capillary bridge breaks fast and leaves 
a well known satellite droplet. By contrast, the plasma 
and the P£^O50i5/85 solution form a viscoelastic thread 
that thins exponentially in time (see Figure [3|. Both so- 
lutions also show the typical blistering instability at the 
very end of the thinning process characterized by the ap- 
pearance of smaller droplets on the filament [18]. These 
results clearly demonstrate that plasma has viscoelastic 
properties. 

Figure [3] shows the minimal width h{t) of the capillary 
bridge as a function of time t — tfinai- Starting at a width 
of approximately hi ~ 0.5 mm, all curves initially de- 
scribe a uniform thinning down to an intermediate level 
of h2 ~ 0.17 mm. At this stage, polymers in the solu- 
tion do not yet affect the flow and the difference between 
the solutions can be only seen during further thinning 
of the filament. While the capillary bridge of the New- 
tonian sample breaks very rapidly (cp. width /13), the 
other samples form a filament which thins exponentially 
over characteristic time scales Ac. Table |T] summarizes 
the corresponding relaxation times as averages over all 
performed measurements. Note that the BSA solution 
does not show any non-Newtonian behavior and there- 
fore we can exclude that the protein surface layer of the 
plasma is responsible for the elastic thinning behavior (in 
agreement with [19j), but some other plasma proteins in 
the bulk that have not yet been identified are responsible 
for the viscoelasticity. This result is also corroborated in 
a cone-plate geometry with a free surface layer: adding 
0.01wt% of the soluble nonionic surfactant (Tween 20, 



critical micellar concentration cmc = 0.07wt% to the 
plasma did indeed suppress the apparent shear thinning 
(data not shown) in agreement with the experiments de- 
scribed in Ref. [5j). Please note that in principle surfac- 
tants can lead to viscoelastic filaments themselves but in 
our case the surfactant dissolved in pure water did not 
show any filament [2QH28] . However, the addition of a 
surfactant to the plasma did not alter the thinning pro- 
cess in the CaBER which again shows that the filament is 
not caused by a protein surface layer. We note that addi- 
tional thinning experiments of the plasma solution were 
performed in which the surrounding air was replaced by 
silicon oil and a pronounced filament was as well observed 
(data not shown). 

Again, we can compare our polymer solutions that 
should serve as model systems with the plasma sam- 
ples. The data presented in Table [T] shows that the 
P£^O50i5/85 solution is a good plasma replacement so- 
lution because it matches the plasma shear and elonga- 
tional properties; by simply increasing the solvent viscos- 
ity of this solution we get the PEObO^^/^^ sample that 
matches the elongational properties of full blood but it 
does not match its shear rheology very well. Therefore, 
we compared a PBS buffer solution with 50vol% hema- 
tocrit without polymer {PBSi^/s^_Hct5o)^ with poly- 
mer (P£^050i5/85-i^ct5o) cind full blood. The polymeric 
hematocrit solutions reflects the elongational properties 
of blood to some extent. However, the PBSi^/s5-Hct50 
solution simply breaks- up like a Newtonian liquid. This 
clearly shows that the non-Newtonian elongational vis- 
cosity (i.e. the elastic properties of the macromolecules) 
of the plasma contributes to the non-Newtonian elonga- 
tional viscosity of blood. Again, an additional clustering 
of the RBC's in the whole blood due to the rouleaux 
formation might explain the remaining discrepancies. 

Next, the effects of this apparent viscoelasticity of 
plasma flow are investigated in a microfluidic contraction- 
expansion device by measuring the pressure drop across 
the contraction part of the geometry [8l|9]. This geom- 
etry has a significant elongation flow component and it 
resembles flow restricted vessels or the flow at a branch- 
ing from a larger to a smaller vessel. The microchan- 
nels are made of PDMS using standard soft lithography 
methods. Their length, width and depth are L = 30 mm, 
W = 400 |LLm and h = 50 \im respectively. These rect- 
angular microchannels exhibit sharp contraction of re- 
duced width w = 25 \im and length / = 100 |Ltm in the 
middle. Note that the ratio of the expansion width to 
the contraction width is 1:16 which sets the (Hencky) 
elongation strain e = In (16) ~ 2.8. Pressure measure- 
ments are performed by directly incorporating two pres- 
sure taps (Honeywell 26PC) into the microchannel, both 
upstream and downstream of the contraction plane. The 
pressure sensors are used to cover a differential pressure 
range of < AP < 56.5 kPa. This type of measurements 
have been successfully applied to investigate the flow of 
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viscoelastic fluids in expansion-contraction microfluidic 
devices 

Results for the pressure drop of both human plasma 
and water as a function of flow rate, or equivalent ly 
Reynolds number (Re) are shown in Figure 4. Note that 
here Re = pVDh/rjo where V is the mean flow velocity, 
p the fluid density, 7^0 the zero shear rate viscosity and 
Dh = 2wh/{w + h) the hydraulic diameter. For exam- 
ple. Figure 4 (inset) shows that the water pressure drop 
increases linearly as the flow rate is increased. On the 
other hand, the blood plasma pressure drop shows a de- 
viation from this linear trend. This nonlinear behavior 
can also be observed by normalizing the pressure drop 
at a given flow rate by the pressure drop of water at the 
same flow rate. Figure 4 shows that the dimensionless 
pressure drop of plasma decreases signiflcantly as the Re 
increases. This indicates an extensional strain-rate thin- 
ning viscosity behavior that has been observed in other 
types of complex fluids such as micellar solutions [29j. 
These measurements clearly show that blood plasma is 
a non-Newtonian fluid. We also checked by flow visual- 
ization that the increase in resistance is not caused by a 
viscoelastic instability, but we found that in principle the 
viscoelastic behavior of plasma could be responsible for 
viscoelastic instabilities, especially in complex flow situ- 
ations and in the presence of red blood cells, i.e. 50% 
hematocrit [30j. 
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Figure 4. (Color online) Dimensionless pressure drop as a 
function of Reynolds number. The data shows higher values 
of pressure drop for human plasma compared to water. Blood 
plasma shows a nonlinear, 'thinning' behavior of pressure drop 
as a function of strain-rate [29 . Inset: Pressure drop as a 
function of flow rate. Blood plasma shows nonlinear behavior. 

In conclusion^ we found a pronounced viscoelastic be- 
havior of human blood plasma in a capillary breakup ex- 
tensional rheometer (CaBER). Similar to dilute polymer 
solution, steady shear flow or small amplitude oscilla- 
tory shear does not indicate any elasticity of the plasma 
because it stretches polymers less efliciently. There are 
several reports on viscoelasticity of plasma in shear flow 



but they all could be attributed to the formation of a 
protein surface layer which is not the case in our mea- 
surements. We could also show that the elongational 
properties of blood are to a large extend determined by 
elongational properties of the plasma proteins. In or- 
der to check the generality of our flndings we prepared 
a plasma replacement fluid with similar elongational and 
shear properties. Differential pressure measurements in 
a microfluidic contraction-expansion geometry conflrmed 
the non-Newtonian character of plasma. Supplementary 
measurements with flow visualization indicate that this 
viscoelasticity might even lead to viscoelastic flow insta- 
bilities. Similarly, it is known that small amounts of 
polymers with elongational properties as our plasma re- 
placement solution can lead to turbulent drag reduction, 
i.e. turbulent instabilities due to high inertia might be 
suppressed by the elongational properties of the plasma. 
Finally, recent numerical results indicate that a slight 
viscoelasticity of the solvent might lead to a pronounced 
cell depletion layer close to the vessel boundaries [31j. In 
view of these flndings, it is expected that the viscoelas- 
ticity of plasma has to be taken into account in future 
studies. 
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